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Abstract 

In this paper, we consider a homogeneous Markov process u) on an ultra- 
metric space Q p , with distribution density f(x,t), x E Q p , t G R + , satisfying the 
ultrametric diffusion equation -^f(x,t) = —D"f{x,t). We construct and examine 
a random variable tz {ui) that has the meaning the first passage times. Also, we 
obtain a formula for the mean number of returns on the interval (0, t] and give its 
asymptotic estimates for large t. 

Introduction 

Ultrametric random processes and their physical and biological applications have recently 
been attracting much attention, especially in connection with modelling the dynamics 
and evolution of complex systems characterized by multidimensional rugged energy land- 
scapes (fitness landscapes) with a huge number of local minima (see, for instance, [T]-[7]). 
It is clear that a description of the dynamics on such landscapes requires adequate ap- 
proximations. As shown recently, a reasonable approximation for the dynamics of some 
biological systems (in particular, proteins) can be chosen in the form of random "jumps" 
between local minima of a landscape, under the assumption that the only key factor is 
the maximal activation barrier on the landscape that separates these local minima [8]. In 
this case, the local minima are clustered in "basins" of minima hierarchically embedded in 
one another. Accordingly, the dynamics of such a system is described in terms of random 
"jumps" between the basins. As shown in recent publications [4] |6] such approximations 
can be naturally described in terms of ultrametric random processes, and it turns out that 
the p-adic pseudodifferential equation of ultrametric diffusion (introduced in [9] and called 
there the equation of Brownian motion on the p-adic line) gives an adequate description 
of protein dynamics [15]. [TP]. 
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Clearly, a physically meaningful application of the ultrametric diffusion equation re- 
quires an answer to some questions pertaining to the description of experimentally observ- 
able quantities. In this connection, it should be mentioned that the observable quantity 
in a real experiment corresponds to some specific degrees of freedom (determined by the 
adopted method of measurement) of a complex system with many degrees of freedom. In 
some cases, the observable quantity corresponds to a characteristic of the system averaged 
over all its states [5], in other situations, it is determined by some specific states [ID]. In 
a situation of this kind, there arise classical probability problems for ultrametric diffusion 
such as the problem of the first passage time distribution and the problem of the number 
of returns. 

In this paper, we consider a homogeneous Markov process £(t;u) on an ultrametric 
space Q p (ultrametric random walk), with distribution density f(x,i), x G Q p , t G R + , 
satisfying the equation 

J^/CM) = -££/(*, t), 

usually called the ultrametric diffusion equation (for definitions and notation see below). 
We consider a specific random process £(t;u)), namely, that whose distribution density 
satisfies the Cauchy problem for the ultrametric diffusion equation with the initial density 
in a domain Z p C Q p . 

Our aim is to construct and examine a random variable Tz p (oj) that has the meaning the 
first time instant when the trajectories of the random process u) return to the domain 
Z p . To study this problem, we first prove that the distribution density of tz p {uj), denoted 
by f(t), satisfies a nonhomogeneous Volterra equation, then we construct a solution of 
that equation and examine its properties. On the other hand, we show that the first 
passage time distribution density can be represented as a functional of a density function 
which is a solution of the ultrametric diffusion equation with the absorbing region Z p . 
It is shown that these two approaches are equivalent. In the last part of the paper, we 
consider the problem of the number of returns to the domain Z p on the time interval 
(0, t\. We obtain a recurrent equation for the probability q^ m \t) of the m-th return on the 
time interval (0, t], as well as a recurrent equation for the probability h^ m \t) of precisely 
m returns on the interval (0,t\. We study the properties of the functions q^ m \t), h^ m \t) 
and obtain a formula for the mean number of returns on the interval (0, t] and give its 
asymptotic estimates for large t. 

Section 1 contains some basic facts from p-adic analysis and the theory of random 
processes. These facts are used for the introduction of the necessary notation and defini- 
tions. In Section 2, we consider the first passage problem for ultrametric random walk. 
In Section 3, we introduce and examine a p-adic analogue of the diffusion equation with 
an absorbing region for the first passage problem. Section 4 is dedicated to the problem 
of the number of returns for ultrametric diffusion. 

1 Elements of ]9-Adic Analysis and the Theory of 
Random Processes 

Let Q be the field of rational numbers and p G Q a fixed prime. Any rational number 
i / can be uniquely represented in the form 
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x = p 1 -, 
b 

where a, b, 7 G Z are integers; a and b are coprime positive integers indivisible by p. The 
p-adic norm \x\ of x G Q is defined by the relations \x\ = p~ 7 , |0| = 0. The completion 
of the field of rational numbers Q with respect to the p-adic norm is denoted by Q p and is 
called the field of p-adic numbers. The set Q p endowed with the metric p(x,y) = \x — y 
is an ultrametric space which is complete, separable, totally disconnected, and locally 
compact. There is a unique (to within a coefficient) Haar measure d p x on Q p which is 
translation-invariant: d p (x + a) = d p x. We normalize this measure by the condition 

dpX — X . 

where Z p = jx G Q p : \x\ p < l| is the ring p-adic integers. There is only one measure 

d p x satisfying the above condition. 

We introduce a class W a (a > 0) of complex-valued functions <p(x) on Q p satisfying 
the following conditions: 



(i) \<p(x)\ < C ^1 + \%\ p \ where C is a constant; 



(ii) there is an integer iV = N((p) > such that ip (x + x') = (p (x) for any x G Q p and 
any x' G Q p such that \x'\ < p~ N . 

A function ip(x) satisfying these two conditions is called locally constant, and N(ip) is 
called its index of locally constancy. For a function (p that additionally depends on a real 
parameter t, we say that (p G W a uniformly in t, if the constant C and the index iV do 
not depend on t . 

Functions in W° with compact support are called test functions (or Bruhat-Schwartz 
functions). The set of all test functions is denoted by D, and the set of distributions on 
D is denoted by D' . 

Let x be a normalized additive character of the field Q p . Then x £ W°. The Fourier 
transform of a function <p (x) G L 1 (Q p , d p x) is defined by 



&{k)= \ x(kx) ip (x) d p x, k G Q p . (1.1) 

Jq v 

For (p (k) G L 1 (Q p , d p k), the inverse Fourier transform is defined by 



(p(x) 



/ x (~kx) (p (k) d p k, x e Q p . (1.2) 

JQr, 



The operator (the Vladimirov pseudodifferential operator |9j), a > 0, is defined on 
functions ip G , < f3 < a, by the formula 



D> (,) = -— -I (1.3) 
1 (-«) ./r)_ a; - yL 



where T p (— a) = 1 ^_ - is the p-adic analogue of the gamma-function. 
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Below, we consider random processes over the field Q p . According to the Kolmogorov 
axioms , a measurable space is a pair {Q, £}, where Q is a set and £ is a a-algebra of 
subsets of Q. A probability space is a triplet {Q, £, P}, where {Q, £} is a measurable space 
and P is a countably additive nonnegative measure on £ such that P(f2) = 1. An element 
A G S is called an event, and the measure P(A) is called the probability of the event A. 
Let {F, B } be a measurable space. A mapping £ : ft — > F is called £|B-meas?ira5/e, 
if £ _1 (B) C E. A E|B-measurable mapping £ is called a random variable with values 
in Y and is denoted by £ = £ (a;). Such a function £ (a>) induces a probability measure 
P$(B) = P{£ _1 (.B)} on sets B G B. The function P^(B) is called the distribution 
function of the random variable £. 

A random mapping of a set T into a measurable space {F, B} is defined as a mapping 
£ (t, uj) : T x ft — > F such that for any fixed t G T it is a measurable mapping from (ft, E) 
to {F, B}, i.e., for any B G B, we have 

{cj G ft : £{t,u) G 5} G E. 

If the parameter t is interpreted as time, a random mapping is called a random process. 

Let Y = Q p , T = B} + . As a probability space one can take Q = Q p and E = B, where 
B is the a- algebra of all measurable subsets of Q p . To define a Markov process on Q p 
homogeneous with respect to time, it suffices to define its distribution density function 
/ (x,t) and the transition density f(y, t\x) satisfying the following conditions: 

1. f(x,t) is B-measurable in x G Q p for any t; 

3. f(y, t\x) > for any x G Q p , y G Q p , and t > 0; 

4. is B x B-measurable in x, y for any t > 0; 

5. Jq f(y,t\x)d p y < 1 for any x e Q p and t > 0; 

6. for any x G Q P , y G Q p , s > 0, and t > 0, the Chapman-Kolmogorov condition 
holds: 

f(z,t + s\x)= f(z,t\y)f(y,s\x)d p y; (1.4) 

7. for any x G Q P , s > 0, and t > 0, the compatibility condition holds: 

f(z,t + s\x)= f f(z,t\y)f(y,s)d p y. (1.5) 

In this case, / (x, t) defines a one-point distribution function for the random process: 

P(B,t)= [ f(x,t)d p x, 

J B 

and f(y,t\x) defines the transition function for the homogeneous Markov process: 



P{B,t\x) = 




( I B (x), t = 0, 
where Ib(x) is the characteristic function of the set B. 
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2 The First Passage Problem 



Consider a homogeneous Markov process £ (t, uj) : R\ x ft — > Q p with the transition 
density 

f(y,t\x) = f(y-x,t)= [ ex P (-\kp) X (k(y-x))d p k. (2.1) 



The function f(y — x,t) satisfies the Markovian conditions 

f(x,t)>0, / f(x,t)d p x = l, 

f(x,t)^8(x) in D' as t -> 0+ , 
f(x-y,t)f(y,t')d p y = f(x,t + t'). 



The function f(y — x,t) of the form ( 12.11) is a fundamental solution of the ultrametric 
diffusion equation 

J^/fot) = -££/(*,*). (2.2) 

This random process was introduced in [9] as a p-adic analogue of random walk (on 
the p-adic line) , and equation (12.21) was interpreted as a p-adic analogue of the diffusion 
equation, although the operator D° is nonlocal and its correspondence to the Laplace op- 
erator is problematic. Note that in contrast to Wiener processes, the p-adic random walk 
£ (t,cj) admits no continuous trajectories, since Q p is a totally disconnected topological 
space. The support of £ (t, uj) belongs to the class of functions without discontinuities of 
the second kind (see, for instance, |9|). The operator D" can be interpreted in more clear 
physical terms, if ( 12 .21 ) is regarded as a kinetic equation [!]-[6], which is justified in view 
of the integral representation (ll.3p of the pseudodifferential operator £)". 

One of the classical problems of random walk on the real line is that of finding the 
distribution function of the random variable describing the first time instant when the 
wandering particle returns to the origin. Consider a similar problem for the p-adic random 
walk £ (£, uj) defined above. 

Let the evolution of the distribution density function <f(x,t) of the random process 
£ (t, lu) be described by the Cauchy problem for the ultrametric diffusion equation 

9 , v 1 f , (p(y,t) - <p(x,t) 

— <p(x,t) = - d p y— — , (2.3) 

with the initial condition 

tp{x,G)=n{\x\ p ), (2.4) 

f 1 \x\ < 1 

where fi(lxL) = < _' p — , ' is the characteristic function of the domain Z v . 

\ 0, \x\ p > 1 p 

Definition. The random variable r Zp (uj) : Q — > R\ defined by the relation 

T Zp (u) = mf{t>0:\t(t,uj)\ p <l, if 31/: \£(1/,u)\>1, < t' < t} 

is called the first passage time of a trajectory of the random process £ (i, u) entering the 
domain Z p (i.e., the first instant when it returns to Z p ). 
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The initial condition (12.41) obviously implies that 

P{u G tt : |f (0,a;)| < 1} = 1- 



Theorem 1. The distribution density function f (t) of the random variable tz p {uj) sat- 
isfies the nonhomogeneous Volterra equation 

9(t)= f g(t-r)f(r)dr + f(t) 
Jo 

with 

Proof. Consider the event A (t^tj) that consists in that a particle staying in the domain 

Q p \Z p goes back to the domain Z p at a time belonging to the interval {ti,tj} (under the 
condition that at t = the particle stays in Z p ): 

A (tj, tj) = 

= {u G Q : St G (ti,tj], hm^(t',uj) £ Z p , Jim Q £(l/,w) G Z p \£(p,u) G Z p }. 

Consider also the event B (U, tj) that consists in that a particle staying in the domain 
Qp\Z p goes back to the domain Z p for the first time at an instant belonging to the interval 
(tj, tj]: 

B (titj) = {lo G Q, : U < r Zp (u) < tj}. 
Let us divide the interval (0, t] into n parts: 

= t < k < t 2 < . . . < tn-l <t n = t. 

We obviously have A{t n - X ,t n ) C U"=i B **)• Since A (*n-i, tn) fl B (tn-i, t n ) = 
B (t n _i, t n ), it follows that 

A (* n _i, t n ) = A (t n _ x , tn) P| ( U 5 U 

\i=l 

n 
i=l 

|J (A (t n -i,t n ) fl £ \ |J ^ (*n-l, ^n) (2.5) 



8=1 



Let P {A (ti-i, ti)} and P {i? (ti-i, tj)} be the probabilities of the events A (t»_i, tj) and 
B(ti-i,ti), respectively. Taking into account (12.50 and the incompatibility of the events 
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B we can write 



P {4 (*»_!,<„)}= (2.6) 



n-1 



^ P {,4 (^_ X , *„) f| 5 (t;_i, **)}+?{£ (*„-!,*„)} = 
i=l 
n-1 

= ^P{A(t n _ 1; t n ) \B (ti^t^P {B (t^ti)} + P {B (t n ^,t n )} = 

i=l 
n-1 

Here, we have used the relation 

P {A (t n _ 1; t n ) \B (U^U) } = P {A (t n _! - U, t n -U)} + e(ti - 

where — — » as — tj_i — > 0. On the other hand, the probability P {A (t n _i,t„)} 
is determined by the solution (p(x,t) of the Cauchy problem for the ultrametric diffusion 
equation (12.31) with the initial condition (12.41) and has the form 

P {A (t n _ ls t n )} = #(t n ) (t n - t n _i) + o (t n - t n _i) , (2.7) 

where g(t) is defined by 



= -T^TT / (2-8) 



r p (-a) Jq p \z p \x\ 

The function is interpreted as the density of the probability to go back to the 
domain Z p at time t. Similarly, the probability P {B (£,_!, ^)} of the first passage to the 
domain Z p on the time interval for the same random process can be represented 

in the form 

P {B (U-i, U)} = f(U) (U - ti-{) + o(U- U-!) , (2.9) 

where f(ti) is the sought density of the probability of the first passage to the domain 
Z p at time U. Now, substituting ( 12.71 ) and ( 12.91) into (12.61 ) and passing to the limit as 
max i=1 n {(ti — tj-i)} — > 0, we obtain a nonhomogeneous Volterra equation of convolu- 
tion type, 

g(t)= [ g(t-r)f(r)dT + f{t). (2.10) 



Note that g(t) is a continuous function, and therefore, equation (12.101) has a unique 
solution in the class of continuous functions (see, for instance, [H]). It is easy to check 
that g(t) is a function with a finite growth exponent for t > 0, and therefore, f(t) has a 
finite growth exponent for t > and there exist Laplace transforms of the functions g(t), 
fit) denoted by G(s), F(s), respectively. Passing to the Laplace transforms in (12.1 Oj) . it 
is easy to find that 
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Let us calculate G(s). Substituting the solution of the Cauchy problem ( I2.3I )- (12.4I) . which 
has the form (see (ED), (Q) 



<p(x,t) = I Q ( | fc | ) exp — \k\pt x(—k x )d p k 



into (I2.8K . integrating the result in x, and then passing to the Laplace transforms in t, we 
get 



B a — 


|A 


la 
\p 


s + | 


k\ 


a 
P 



G(s) = I n(\k\ p ) / 1 J d p k = (B a + s)J(s)-l, (2.12) 

where 



[i-p- 1 : 
i — p~ 



B a = : JiJv P~ a <B*<l, 



oo 

J (a ) = (l-p-i)5>-__. (2.13) 

n=0 6 + P 

Substituting (12 . 1 2H into (12 . 11 h . we obtain the Laplace transform of the desired function: 

F(s) = 1 - 77^ \ w x • (2.14) 

The function F(s) is analytic in the domain Res > and tends to zero as \s\ — > oo, 
uniformly with respect to arg s. The function F(s) is the Laplace transform of the function 
f(t) with zero growth exponent: \f(t)\ < M. Now, it is not difficult to show that fit) 
has the following properties: 

1. For a > 1, we have J °° f(t)dt — F(0) = 1, which means that for a > 1 the particle 
is sure to return to the initial region, and therefore, on an infinite time interval will go 
back to that region infinitely many times. In this case, however, there is no finite mean 
waiting time for the first passage: 

(r Zp ) = tf it) dt = - lim -F(s) - +oo . 



Res>0 



2. For < a < 1, we have / °° f(t)dt = F(0) = ^ = C Q < 1. This means 

that for small a, there exist trajectories of the unltrametric random walk that abandon 
the initial region never to go back. Note that for the real-valued Brownian motion the 
return property of its trajectories is missing only if the dimension of the space is greater 
than two. 

Consider more closely the function F(s). Clearly, it has simple poles at s = —\k, 
k — 0, 1, 2, . . ., which are simple roots of the equation J(s) = 0, and s = —B a = — A_i. 
From (12.131) . it is easy to see that the values belong to the interval p~ a ( k+1 *> < \ k < p~ ak . 
The point s = is essentially singular and is a limit point of the poles. The function F(s) 
is non-meromorphic on the complex plane, and this is an obstacle to the application of 
the residue theory for the calculation of the inverse Laplace transforms. To overcome this 
obstacle, we first prove the following result. 
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Lemma 1. The function F(s) can be represented as an infinite sum of terms with simple 
poles at the points s = —Xk, k = — 1, 0, 1, 2, . . namely, 



m = E 



k=-l 



S + X k 



(2.15) 



where b k are the residues of F(s) at the points —Xk- On any closed set G that does not 
contain s = 0, the series (12.151) becomes uniformly convergent, if its finitely many terms 
with poles in G are dropped. 

Proof. Consider the auxiliary function 

1 



w 



lim $(0) = lim F(s) = 0, lim $(w) = F(0). 



The function $(«;) is analytic on the complex plane except at the simple poles = — 
w k = —j- , k = 0,1,2 — By the Mittag-Leffler theorem (see, for instance, [32]), <3>(u>) 
can be represented in the form 



$(w) 



E 

fc=-i 



Ck 



w -w k 



Pk + c, 



where c is a constant and c k are the residues of and this series becomes uniformly 

convergent on any closed bounded set, if its terms with poles in that set are dropped. 
Then 



F(s) 



k=-l \ s A fc / k=-l v k 

oo 

E 



Pk 



k=-l 

and since $(0) = 0, we have 



— — — + CfcAfc - p k + c, 
s + Xk 



00 \ 2 

F{s) = E CkK 



k=-l 



s + X h 



Letting c k X\ = b k , we obtain (12.151) . 

From Lemma 1 and the uniform convergence of (12.151) in the domain Res > Sq > 0, 
we see that for the calculation of the original function f(t) it suffices to apply the inverse 
Laplace transformation to the series (12 . 151) term-by-term. Thus, we get 



f(t) = L;i t [F(m) = £m^ 



fc=-l 



s - Xi 



and finally, 



f(t) = £ b k exp[-X k t], 



(2-16) 



k=-l 
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where 

6 -> - jRsy (2 - 17) 

6 fc = _^ fc = 0,1,2,.... (2.18) 

It is not difficult to see that the series YlT=o^ k * s convergent and majorizes the series 
(I2.16p . which implies uniform convergence of the latter and the continuity of f(t). 
The above results can be summed up as follows: 

Theorem 2. The distribution density for the first passage times of a trajectory of the 
ultrametric random walk can be represented as a uniformly convergent series (12.161) whose 
coefficients are defined by (12.171 ) and ( 12.181 ). 

Let us go on with the examination of f(t). It is not difficult to show that 

1) lim /(f) = lim sF(s) = 0; 

2) lim f(t) = lim sF(s) = 0. 

Res>0 

Then, since the function f(t) is positive and continuous, it must have a maximum. 
Let us show that this maximum is unique. 

From the first limit, we have YlT=-i = 0. Thus, the series ( 12.161 ) can be represented 

as the difference of two monotonically decreasing strictly concave down functions, fit) = 

~~/i(f) + /2(f) > 0, and therefore, sup |/i(f) — f2(t)\ is unique. 

teR+ 

The asymptotic behavior of the function f(t) for all a is described by the following 
theorem. 

Theorema 3. For the first passage time distribution density fit) the following estimates 
hold: 

A(a)t-^ (1 + o(l)) < f(t) < B{a)r^ (1 + o(l)) for a > 1; (2.19) 
A(a)r« (1 + o(l)) < f(t) < B(a)t~ « (1 + o(l)) for a < 1; (2.20) 

A(a)-^- (l + o(l)) < f( t )<B(a)-^- (l + o(l)) for a = 1, (2.21) 

where o(l) —> as t —>■ oo, and A(a), B(a) are functions of a and p. 
This theorem is proved in Appendix B. 



3 p-Adic Analogue of the Diffusion Equation with Ab- 
sorbing Region for the First Passage Problem 

For the classical problem of random walk of a particle on a straight line, it is well-known 
that the distribution density function for the first instant at which the particle reaches a 
given domain can be found from the solution of the diffusion equation with an absorbing 
region (see, for instance, [H]). We are going to show that a similar approach can be 
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used in the p-adic case: the distribution density function f(t) for the time of the first 
return to the domain Z p can be obtained from the solution of the Cauchy problem for the 
ultrametric diffusion equation with the absorbing region Z p , i.e., the equation 



di/j(x,t) 1 / f ip(y,t) - ip(x,t) . . f ip(y,t) 




dt r p (-a)\Jn \x-y\p 1 



d P y-n(\x\ p ) ! ^ {v V a+1 d P y), (3.1) 



with the initial condition ij)(x, 0) = fi(|x| p ). The second term in the right-hand side 
of equation (13.11) is equal to the probability of transition from the region Q p \Z p to the 
absorbing region Z p per unit time. Since this transition for all trajectories of the random 
walk (|3.1I) is always the first one, it follows that the probability density of this passage at 
the instant t is defined by the formula 



r p (-«) J Qp \ Zp \x\ 



/(*) = ~ r F~T~~\ I ^dx. (3.2) 



Thus, we have two approaches to finding a solution of the first passage problem. Their 
equivalence is established by the following theorem. 

TeopeMa 4. The first passage time distribution density function obtained from the solu- 
tion of the Cauchy problem for the ultrametric diffusion equation with the absorbing region 
Z p coincides with the solution of the Volterra equation (12.101) 

Proof: Let us apply the Fourier transformation to ip{x,t) with respect to the p-adic 
variable x and then the Laplace transformation with respect to the real variable t. Denote 
the resulting Fourier-Laplace transform by ^ (k, s). From (13 - 1 ft . taking into account the 
initial condition i/j(x, 0) = f2(|x| p ), we obtain the following nonhomogeneous Fredholm 
equation for \I> (k, s): 

8*(k, s) = Q(\k\ p ) - \k\ft(k, s) - Q(\k\ p ) [ *{q, s) (B a - Q(\q\ p )d p q, 

JQp 

or 

*(*, s) = _ !S| f ^ s) {Ba _ K) n(\q\ p ) dp q. (3.3) 

S + \K\ p S + \K\ p J Qp 

Multiplying equation (13. 3h by ( B a — \k\ p ) VI ( \ k\ ) and integrating the result, we get 



! *(k, s )(B a -\k\;)n(\k\ p )d p k= f n(\k\ p ) - ,'5 ^ + 

f B — \k\ a f 

+ / "(W \li\ a P d P k / ^(q,s){B a -\q\ a p )n(\q\ p )d p q. (3.4) 
JQp s + \ K \p JQp 

Note that J Q \I/(g, s) (£> Q — \q\ p ) fl(\q\ p )d p q = F(s) is the Laplace transform of the first 
passage time distribution density function f(t) defined by ( 13.21) . Now, in view of (12 . 1 2f) . 
we can rewrite equation (13.40 in the form 

F{s) = G(s) +G(s)F{s). 

Comparing this with (|2.11l) . we see that the solution of the last equation coincides with 
that of the Volterra equation (12.1 Oj) . 
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4 Number of Returns for Ultrametric Diffusion 



In this section, we consider some questions pertaining to the probability of the m-th return 
on the time interval (0, t) and the growth of the number of returns with the growth of t. 
For the probability space {ft, E, P}, consider a random process 

N Zp (t, uj) : ft x N C Z+, N Zp (0, uj) = 

that describes the number of returns of a particle to the domain Z p on a finite time interval 
(0, t] . Let us calculate the probability of the m-th return of a particle to Z p on the interval 
{0,t}. Consider the event Q™ = {uj G ft : N Zp {t,uj) > m} that consists in that a particle 
staying in the domain Q p \Z p goes back to Z p for the m-th time at an instant from the 
interval (0, t], or equivalently, that a particle visits the domain Z p at least m times on the 
time interval (0, t]. Denote the probability of this event by P{Q™} = q^ m \t). Obviously, 
= for all m > and Q° t = ft for all t > 0. 

TeopeMa 5. The probability q( m '(t) of the m-th return on the interval (0, £] satisfies the 
recurrent equation 

q {m \t) = j\ {m - l \t-r)f{r)dr, m > 1 , 
g(°>(t) = l, m = 0, 
where f(t) is the distribution density for the first return time. 

Proof. This statement is proved along the same lines as Theorem 1, and therefore, we 
just outline the main steps. 

Consider the event Q™. Let B T+dr = {uj 6 ft : r < r Zp (uj) < r+dr} be the event of the 
first return to the domain Z p on the time interval (r, r + dr}. Then, Q™ C U r e(o t] B T +dr 
and for m 7^ we can write 

q? = q?C\ \ U B ^r) = U (^rD 5 ^ 

\Te(o,i] J r£(0,t] 

Since the events B T+dr are incompatible for all r G (0, i], we have 

re(0,t] r6(0,t] 

Observing that P{Q™ \B T+dr } = PjQ^ 1 }, we obtain 

HQ?) = E pi^+^jpigr. 1 }, m > i, 

re(o,t] (4.2) 
P{Q° t } = 1, m = 0. 

Finally, recalling that P{Q?} = q {m \t) and using the symbolic formula P{B T+dr } = 
f{r)dr (its meaning is clear from the rigorous arguments in the proof of Theorem 1), we 
obtain the desired recurrent relation (14.11) . 
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For the Laplace transforms, equation (14.11 ) reads 

QM(s) =Q( m - 1 )(s)F(s) 



Q (0) (s) = s- 



i 



(4.3) 



where Q ( s ) is the Laplace transform of q^ m \t). 
The solution of ( 14.31) has the form 

Q m (s) = - s (F(s)) m , (4.4) 

where F(s) is defined by (12. 141) . 

Using the recurrent equation (14.11) and the properties of / (t) and F(s), it is not 
difficult to show that the functions q( m \i) have the following properties: 

1. Each g(™)(t) is a monotonically increasing function oft; the function j^q^it) > 
has the meaning of probability density for the m-th return. 

2. Iim4-g (m) (t) = 0, lim ^q {m \t) = 0. 
t-+o dt y 1 t-*oo dt v 1 

3. limg( m )(t) = 0. 

4. lim g( m )(t) = j « < l' Reca11 that the <l uantit y G « = ■§* (^t) < 1 

(see Section 2, property 2 of /(£)) is the measure of return trajectories of the ultrametric 
random walk, which is the same as the probability of the first return on the infinite time 
interval. 

5. q( m \t) < g^ m ~ 1 ^(t), m > 1, i.e., the sequence q^ m '{t) is monotonically increasing 
with respect to m for any t. 

With the help of the recurrent equation 1 14. 1R and properties 1, 2 of the function q^ m \t), 
it is not difficult to show that the return probability densities have a maximum, which is 
unique, and thus, we have a single-mode distribution. 

In view of property 4, it is only for a > 1 that one can speak about the mean waiting 
time of the m-th return. For a > 1, the m-th (in particular, the first) return is a certain 
event, but its mean waiting time is infinite. Indeed, 

r^q {m \t)dt= lim (-m(F(s)) m_1 ^F(s) J =oo, 



since lim 4-F(s) = — oo. 



+0,Res>0 ds 



Next, consider the problem of finding the probability of precisely m returns on the 
time interval (0, t]. Let H™ = {u 6 fl : N Zp (t,u) = m} be the event that on the time 
interval (0, t], the particle goes back to the region Z p precisely m-times. We are interested 
in the probability of this event, P{H™} = /i (m) (t). 

TeopeMa 6. The probability hS m \t) of precisely m returns on the time interval (0, t] 
satisfies the following recurrent equations : 



hP\t) = 1- / f(r)dr, m = 0. 
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where f(t) is the probability density for the first return time. 
Proof: We obviously have H^ 1, = QT\QT +1 , an d therefore, 

h^ m \t) = q {m \t) - q {m+1 \t). (4.6) 

Substituting f l4.ll) into f 14.6l> . we obtain the recurrent equations ( 14.51) . The theorem is 
proved. 

Let us examine more closely the probability distribution function for precisely m re- 
turns. In terms of Laplace transforms, equation (14.51) has the form 

ifM(s) = H ( - m - 1 Hs)F(s), 

(4.7) 

#(°)( s ) = s- 1 (l-F(s)), V ' 

where H^ m '(s) is the Laplace transform of h^ m \t). From (14.71) . we obtain the following 
expression for the Laplace transform of the solution of equation (14.51) : 



H (m \s) = -{l-F{s)){F{s)) m . (U 
s 



Using (12.141) and the solution (14.81) . it is not difficult to establish the following properties 
of^ m )(t): 

1. h^ m \t) is a positive function such that 

lim h {m) {t) = lim (1 - F{s)) {F{s)) m = 0, a > 1; 

t^O s^oo 

2. \imh( m \t) = 0, a<l, lim h (m \t) = (1 - C a ) (C a ) m , a < 1 ; 

3. /i (m) (t) has a maximum, which is unique; 

4. h {m \t) < U m - l \t). 

What is the mean number of returns fi(t) on the time interval (0, i]? Usually, it is 
expected that the mean number of returns is proportional to the walk time. By definition, 
we have 

oo 

/i(t) = X><">(f). (4-9) 

n=l 

Theorem 7. The mean number of returns on the time interval (0, t] is determined by 
the formula 

fi(t) = [ g(r)dr, (4.10) 



o 



where g{t) is defined by (12.81) and is the density of the probability to return to the domain 
Z p at the instant t. 

Proof. Writing the expression (14.91) for Laplace transforms and using (14.81) , we obtain 

oo 1 oo 

M( s ) = j2 nH(n \ s ) = -( i - F (s))Y, n ( F ( s )) n - ( 4 - n ) 



n=l n=l 
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Since |-F(s)| < 1 for Res > 0, the series in ( 14.111 ) can be summed and we have 

M(s > = 1T%> ' < 4 - 12 ' 

Hence, using (I2.1ip . we get 

M(s) = -G(s). (4.13) 

Applying the inverse Laplace transformation, we obtain (14.101) . 

Now, let us calculate the average number of returns on the time interval (0,t], using 
(14. 1 Oh . Integrating equation (12. 3ft over the domain Z p and taking into account (12.81) . we 
obtain the following equation: 

l Szp (t) = -B a S Zp (t)+g(t), (4.14) 

where Sz p (t) = f z <p(x,t)dx. From the solution of the Cauchy problem for equation (12.31) 
and the initial condition (12.41) . we have the following expression for S Zp (t): 

S Zp (t) = (l " -) f>~ n exp [-p- an t] . (4.15) 

V P ' n=0 

The series (14.151) is uniformly convergent, and therefore, using (14.14p . (14.151) in (14.10p . we 
easily obtain 

/*(*) = B a (\- -) f^P^ (1 - exp [-p- an t]) + 

^ P' n=0 

+ f 1 ~ ~) £p~" ex P [-P~ an t] - !• (4-16) 

^ P' n=0 

Note that the first series in ( 14.161) is convergent for all a > 0, although there is 
no uniform convergence for a > 1. The asymptotic behavior of the second series is 
characterized by the function t~& (see formula (1A.21) in Apendix A). 

Theorem 8. The following asymptotic estimates hold for the function fji(t): 

V P J \a J [a — I) mp 

<fi(t) <p(l--)r(-) Ba t^+0(ri), a>l, (4.17) 
V P J \a J (a — 1) mp 



V 



1 (l--)^ L T{l)lnt + 0{t- 1 ) < fM(t) < 
\ p J mp 

< p ( 1 — - J -^r(l)lnt + 0(t- 1 ), a = l, (4.18) 
\ p J In p 



i-c a 

1 — C a * J V p7 \ a J alnp 



< T %- " P 1- " fl - zV f —1 7^r*-^(l + or 1 )), (4.19) 

< a < 1, 
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where T(x) is the gamma- function and C a = ^ {^p-. T ) ^ s ^ e probability of the first 
return on the infinite time interval. 

To obtain the asymptotic estimates of fi(t) for a > 1, it suffices to integrate the 
asymptotic estimates 

P -\1 - p- X )B a V -^-r-« < git) < p(l - p- x )B a T -^±t^ 
amp amp 

obtained from an asymptotic estimate for the series S(t) (see formula ( 1A.2I) in Ap- 
pendix A). This integration is justified, since the function g(t) and its asymptotic bounds 
continuously depend on t > a, a > 1, are strictly positive for large t, and 

oo ^ 

t~^dt = +00. 



a 



To obtain the asymptotic estimates of for < a < 1, it suffices to rewrite (14.160 
in the form 

= Bjl-^jTp- (1 - a)n (l-e W [-p- an t}) + 

+ jTp~ n exp [-p- an t] - 1 = 

V P J n—n 



n=0 

00 



\ P' n =0 

and use formula ( 1A.2I) from Appendix A. 

Appendix A 

Here, we obtain an asymptotic estimate for the series 

00 1 

S(t) = T a ~ ie ~ (6) ~ U ' t>0, keN, a>l, b>l, (A.l) 
i=o z 

for £ > 1. 

Lemma Al. For £/ie series (I A. II ) £/ie following estimate holds for t 1 : 
(lno)^ 1 (ln(6t))~ fc (6t) - ^ T f ^ ] (1 + o(t)) < S(t) < 



In 6 



In a 



< a (In o)*" 1 (ln(t))- fe (t)HS r ( ^ ) (1 + o(t)) , (A.2) 



where T (z) is the gamma-function. 
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Proof. Note that -^a x is a decreasing function and e ^ 3=4 is an increasing function of 
x. Therefore, on the interval « < i < i + 1 we have the inequality 

±-a- x e- b - {x - lh < a- l e- b ~ H < -J—a^e"*-'*. (A.3) 
x" (x — If 

Integrating (1A.3I ) in x from i to i + 1, we get 

■i+l i . /-i+l 



a 



1 / ±a-t x -%- bH '- lh dx < a-'e-^t < a / -. a^e^dx. (A.4) 

A Ji (x - If 

Now, summing the inequalities (IA.4I ) with respect to i from to oo, we find that 

f°° 1 

S min (t) = a' 1 I - a ~ (X ~ 1)e ~ b ~ (X ~ lhdx ^ S ^ ^ 



x 



<a I - a x e dx = S mSLX (t), 



where 



k-l n_/t.\\-k / t i \ t~t 



S^t) = (Inb)^ 1 (\n(bt)) 



S max (t) = a (In bf- 1 (ln(t))- fc jf (l - "* yfef-^dy. 



y 2 ~ 1 e-^ = 7 (fc) (^,^) 



Let 

In x , 

and note that lim x k > (z,x) = f?° t z ~ 1 e~ y dy = T(z) and (z,x) = \^.t z ~~ x e~ v dy 

x^oo 

j(z,x) . Then, for x 1, we can write 



S mhi {t) = (mfe)"- 1 (Hbt))~ k (bty 1 ^ r I pi 1 (1 + o(l)) 



In a 
mT 



S max (£) = a (ln&f- 1 (ln(t))-* T ( ^ ) (1 + o(l)) 



Info 



and therefore, the estimate ( 1A.2II holds. 



Appendix B 

Proof of Theorem 3 from Section 2 

To estimate the function /(t), we first estimate the coefficients of the series (12.161) . 
These coefficients coincide with the residues of the function F(s) at the poles s = — A&, 
fc = -1,0,1,2,... (see (I2T7D . (BHD ): 

, 1 r 8 + A fc _ 1 1 1 

fc " (S a - A fc ) J(s) (Sa-A fc ) J'(-Afc) (5 a -A fc ) UA; ' 
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where 



it*. 



-i 



oo 

E 



P~ 



-1 



(B.l) 



Recall (see Section 2) that the poles s = —X k , k = 0,1,2,... coincide with the simple 
roots of the equation Y^=o P~ n s+ p- a n = and the values lie on the interval p- a ( k+1 ) < 
X k < p~ ak ■ The point s = is a limit point for the set of poles. Let us examine the 
behavior of the poles and the residues Uk for large k. We pass from A^ to new variables 
8 k , setting 

A fc = p~ ak (p~ a + (1 - p- a )5 k ) = p- Q ( fc+1 ) + p- ak (l - p' a )5 k , < 4 < 1 
and let 

p-r> 

(m—aik+i) _|_ n—ak 

1=0 



v k = (l -p *) u k 



It can be shown that flB.31) implies the following inequalities for z/ fc : 



-an\2 



-1 



(B.2) 



(B.3) 



1 _ p (2a-l) 



p (2a-l)* [ (j9 -a + (1 _ p -a )4) 2 (1 _ p2 a-l ) + (1 - p""^ 

_9 n -1 



P 



+ 



(B.4) 



P" 



1 



<V k< (2f 



+ (1 - p- Q ) 2 (p"« + (1 - p-")^) 2 (1 " P- 1 ) 

Since A&, = 0, 1, . . ., are zeroes of the function Y^=o s -p- an ' we have 



■p(i- p -") 2 *|. 



E 

71=0 



P" 



P" 



*(*+!) + p-o*(l - p- a )5 k ~ P~ 



o. 



which implies the following estimate for 8 k : 

Ap 2 ~ a a k 



1 p + 1 + p 2 a a k 



a k 2p 



,2-a 



1/2 



(p+l+p 2 - a a k f 



< 



<5 k < 



1 l-p- a - 1 

V k (p-i){i-p- a ) 



(B.5) 



where ak = k for a = 1 and a k = 1 ti^ 4 " } for a ^ 1. The quantities have the 
following asymptotic behavior for — > oo: 



a fc = - — + o(l) for a > 1, 
1 p 



pi 

l >l " ^-^(1 + 0(1)) for a< 1, 



«fc = — 7P^ 
a — 1 

ak = k for a = 1, where o(l) — > as — > oo. 
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Using these relations and ( 1B.4I ). ( 1B.5I) . it is not difficult to obtain the estimates 
D{a) (1 + o(l)) < 4 < 17(a) (1 + o(l)) , 

D{a)p {l ' 2a)k (1 + o{\)) <u k < U(a)p {1 - 2 ^ k (1 + o(l)) for a > 1, (B.6) 
D{a)p (a ~^ k (1 + o(l)) < 5 fc < f/(a)p (Q - 1)fc (1 + o(l)) , 

D(a)p~ k (1 + o(l)) < M fe < {7(a)?T fc (1 + o(l)) for a < 1, (B.7) 
D^AT 1 (1 + o(l)) < 4 < {/(a)AT 1 (1 + o(l)) , 

b{a)p- k k~ 2 (1 + o(l)) < M fc < U(a)p- k k~ 2 (1 + o(l)) for a = 1, (B.8) 

where -D(a), {7(a), 7) (a), {7(a) are functions of a and p whose expressions are too lengthy 
to be written out here. 

Let fit) be the probability density function ( 12.1 6[) for the first passage times. Taking 
into account the above notation, we can write 



f(t)= b k ex P (— A**) = — exp (-£ a t) + 

fc=-i Jl aJ 



+ r =^T^ ex P + -P" a ) 5 o)i] + <?(*), 

B a -P 01 -{I- P a )0 

where we have set 

»W s ^. r ^. r11 .^ -' [-Or™ + „-(! - „-)&)«] 

For #(£), we have the estimate 

OO , oo 

— J] Mfc exp [-j9^ fe t] < <#) < - - ^^exp • 

Q fe=i a ^ fe=i 

Further, taking into account (1B.6I) - (1B.8I) . we find that 

EQ\£ p v-*>.)k (i + o( i)) exp < ^ (t) < 



a k=l 



< U ^_ a P [1 ~ 2a)k (1 + o(l)) exp [-p-Vt] , a > 1; (B.9) 
a — ^ fc=i 



^f>- fc (l + o(l)) exp [-p-*t] <<?(*)< 



5 " 

< a Ep' fc (l + o(l))exp[-p-V a t], a<l; (B.10) 

^ a p fc=i 

E tt ( x + °W) ex P [-p~**] < < 
1 k=i 

1 p k=i 

19 



From (jR9j>- (jBHl) . using the inequalities (TOl ). ([Oil , we obtain ( 12191 ) - (I2T2B . The proof 
is complete. 
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